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We develop theoretical and numerical tools for the quantification of entanglement in systems with 
continuous degrees of freedom. Continuous variable entanglement swapping is introduced and based 
on this idea we develop methods of entanglement purification for continuous variable systems. The 
success of these entanglement purification methods is then assessed using these tools. 



Pacs No: 03.67.Hk, 42.50.-p 



INTRODUCTION 



Recently the theoretical idea of teleportation in sys- 
tems with continuous variables has been developed Q-^] 
and shortly afterwards demonstrated by Furusawa et. al. 
H . The efficiency of entanglement manipulation proto- 
cols critically depends on the quality of the entanglement 
that one can generate. It is therefore essential firstly, to 
be able to quantify the amount of entanglement in sys- 
tems with continuous variables and secondly, to develop 
methods of entanglement purification that allow the dis- 
tillation of entanglement by local means in such systems. 
The present paper provides answers to both of these es- 
sentials. 

Most analysis of entanglement in continuous variable 
systems relies on expressing the state of the system in 
terms of some discrete but infinite basis, often with math- 
ematical techniques from quantum optics y,Ql- In this 
way the above procedures can be demonstrated and the 
entanglement quantified. Such analysis is convenient but 
not necessary for the theoretical description of these pro- 
cesses. However, for quantifying entanglement it is es- 
sential. 

In pure finite bipartite systems the Schmidt basis of a 
particular system is useful, particularly because its coef- 
ficients allow us to calculate the entropy of entanglement 
H. In this paper we will first describe how such a ba- 
sis occurs in continuous systems from the mathematical 
area of integral eigenvalue equations and how the entan- 
glement can be calculated from it. We then present two 
classes of continuous entangled states and calculate their 
entanglement using a numerical solution from statistical 
mechanics g]. 

Purification is the process by which the entanglement 
in a bipartite state shared between two spatially sepa- 
rated parties (Alice and Bob) can be increased using only 
local operations and classical communication \^. There 
are many methods of achieving this for discrete pure sys- 
tems such as the Procrustean method [Q (for one copy 
of the bipartite state) and using entanglement swapping 
| [10{ (using two copies of the bipartite state one of which 
is not shared). So far it has not been shown that any 
of these procedures have a continuous analogue which is 



able to purify an entangled state. 

Here we will show that a continuous generalization of 
an entanglement swapping procedure using the continu- 
ous controUed-NOT and Hadamard gates introduced by 
Braunstein pTl is able to produce purification in one of 
our two classes of continuous entangled states. 

Section | first presents necessary mathematics from 
integral eigenvalue equations in the context of contin- 
uous quantum mechanical systems. Section O gives two 
quantum gates and presents the two classes of entangled 
states. The calculation of the entanglement is attempted 
analytically as for as possible in section [II and the nu- 



merical procedure is presented and used where necessary. 
In section [V entanglement purification is attempted and 
validated using the techniques developed in the previous 
section, and finally conclusions are given in section M. 



I. PRELIMINARY MATHEMATICS 

In this section we present some of the mathematical 
tools that are used in this paper for the description of 
continuous variable systems. For continuous states we 
can express a pure bipartite system as 



ip(x,y)\x) 



dxdy, 



(1) 



where \x) are the eigenstates of the continuous system 
(position, say). We wish to find the Schmidt basis for 
the bipartite system for which either partial density op- 
erator of the system, e.g. 



Pi = / 2(2; 1-0) 12 i2(V'|a:)2 dx 
pi{x,y) \x) 11 {y\ dxdy 



(2) 



is diagonal. All the necessary mathematics is covered in 
the area of integral eigenvalue equations ||l^] : so suppose 
that we wish to find the eigenvalues of the kernel pi , that 
is, we wish to find 0i(a;)'s such that 



pi{x,y)(j>i{y)dy = A.,0i(a;), 



(3) 



where A^ is the eigenvalue corresponding to (j)i{x). For 
Hermitian kernels (for which pj(y,a;) = pi{x,y)) the 
eigenvalues are real and the set of eigenfunctions will be 
linearly independent, complete and also orthogonal 



i{x)(i)*{x)dx = 0, i T^ j. 



(4) 



A special class of kernels are those which are quadrat- 
ically integrahle: 



K{x, y)dx dy converges, 



(5) 



where the ranges may be finite or infinite. These are 
called L2 kernels as they are integrable over the L2 space. 
Some basic properties of Hermitian L2 kernels are: 

• they have infinitely many non-zero eigenvalues or 
are PG (Pincherle-Goursat) kernels, ones which can 
be decomposed into the form 



P{x,y) ^^X,{x)Y^{y), 



(6) 



where {Xi{x)} and {Yi{y)} are two linearly inde- 
pendent sets of functions. 

• Their eigenvalues have no accumulation points, i.e. 
the eigenvalues do not form a continuous set, ex- 
cept at zero. 

• The set of eigenvalues converges in the following 

ways 

^ Ar = A„ = Trip",) = J Pi (^, ^) dx 

V finite n, (7) 



2 = 1 



where 



Pii^-.y) = / pi{x,z)pi{z,y)dz 



(8) 



and further powers of pi can be obtained by induc- 
tion. 

• A Hermitian L2 kernel can be written as the ex- 
pansion 



Piix,y) =^A,(/),;(x)(/)*(y) 



(9) 



(where the (j)i (cc) are normalised to the square mod- 
ulus) or the kernel is PG when there are only n 
terms, with n as in equation (ph. 



We see from the above that for Hermitian L2 kernels 
we can always find a diagonal Schmidt decomposition 
into an orthogonal basis given by the kernel's eigenfunc- 
tions. The dimension of this basis therefore depends on 
the number of linearly independent eigenfunctions the 
kernel possesses. 

The measure of entanglement for a pure bipartite sys- 
tem is now easily generalized to continuous variables and 
is just the Von Neumann entropy of either partial density 
operator of the system 



E{pi2) = S{pi) = 5(^2) = - ^ \il0g2\i, 



(10) 



the number of terms in the sum depending on the form 
of pi or p2. We will call this the entropy of entanglement 
or just the entanglement. 

Properties such as concavity, subadditivity, strong sub- 
additivity and the triangle inequality also follow for this 
measure of entanglement as they do its finite partner 
pa provided the relevant quantities converge when we 
deal with infinite systems. The invariance under unitary 
transform of the subsystems can also be proved, where 
the transform is of the form 



U\x) 



U{y,x) \y) dy 



(11) 



with 



C/tC/|x) = |a;) ^ U*{y,z)U{y,x)dy ^5{x-z). (12) 



Transforming the subsystems independently 

p'i2 = U1U2P12UM (13) 

and tracing out system 2 gives 

p[ = Umul (14) 

as the trace is invariant under unitary transform. It can 
then be shown that p', has the same eigenvalues as pi 



/ 



pi {(pi), = Xi {(pi), or 

pi{x,y)4>,{y)dy ^ X,(t)^{x) (15) 



but its eigenfunctions are related to those of p, by 

U^)^ fu^{y,x)^',{y)dy. (16) 



These mathematical techniques, particularly the latter 
method of showing the equivalence of eigenvalue equa- 
tions by transforming the eigenfunctions, will be used in 
the next sections to determine the amount entanglement 
in continuous variable systems. 



II. QUANTUM GATES AND CLASSES OF 
ENTANGLED STATES 



We now turn to some classes of entangled states in 
continuous systems, which will be a convenient point to 
introduce two quantum gates generalized from discrete 
gates [O. The first is the continuous Hadamard trans- 
form which is in fact just the Fourier transform, and in 
which we will include the scale length, a^ explicitly: 



T\x) = 



1 



Jlixyja 



\y) dy 



(17) 



This is also, of course, the transform used to go from the 
position to the momentum basis if we set a = 1 and work 
in units Ti = 1/2. The inverse, T^ is obtained by replacing 
i by —i giving the result that TT^ \x) = JF^JF |a;) ~ \x). 
Note that the scale length, a, is normally inserted to 
make the expression in the exponential dimensionless but 
we will include it as a convenient scale with which to com- 
pare various lengths in the states we are about to form. 
The second important gate is the two particle 
controUed-NOT gate: 



Ci2|2:)il2;)2 = \x)i \y + x) 



2 ' 



(18) 



which is not its own inverse. This is obtained by replac- 
ing the -|- with a — sign on the right hand side. 

Another version of the controUed-NOT, which is its 
own inverse, is 



Such states can be used to demonstrate teleportation 
of an unknown state, the fidelity of the teleportation in- 
creasing as a and /3 — > 0, where the state becomes like an 
infinitely squeezed two mode squeezed state or an EPR 
state Um. We will call the states \Bai3{xi,X2))i2 partially 
correlated entangled states. 

The second class of states we will call two-mode cat 
states |l^,|l^, as they are states which are like two 
Schrodinger cat states whose locations are correlated 
with each other quantum mechanically: 



\Cid)), 



Aoe 



-{x~df~{v+df 



Aie 



-(x+df^(y-df 



\x)i \y)2 dxdy 



J=0 



j^.^~'{x-(-iydY-{y+(--i-ydY 



\x)i \y)2 dxdy. 



(23) 



The complex coefficients, Aj, are such that |Aop -I- 
|Aip = 1 (so the state is not normalised correctly). This 
state is a superposition of the first particle being located 
around d and the second around ~d and vice versa and 
so is not of great use in teleportation. The scale length 
does not appear here (it is set to unity) as an increase in 
scale length is equivalent to a decrease in the value of d. 
It is now natural to ask what the amount of entangle- 
ment in the states \Ba0{xi,X2))i2 and \C{d))■^^2 is, which 
is the aim of the next section. 



C'\x,y) 



X + y X ~ y 



V2 ' V2 



(19) 



III. QUANTIFICATION OF THE 

ENTANGLEMENT 



This is a more symmetric gate and is in fact the trans- 
form produced by a 50:50 beam splitter on the quadra- 
ture wavefunctions in quantum optics M. However, we 
will proceed with the original definition of the CNOT for 
simplicity (primarily of mathematics). 

We can now define the 'entangling' operation and its 
inverse: 



'12 



■ CuJ^i 



'-'12 ^ ■^l'-'12' 



(20) 



and form our first class of entangled states by applying £ 
to two Gaussian wavepackets (aside from normalisation) : 



{x- Xi)^ 



\Ga{xi))^ = / exp 
and \Gi3{x2))2j resulting in the state 

Cl2.Fi|G„(.Ti))JG^(x2))2 
" 1 / 

exp 



dx 



(21) 



2 2 

~x a 



1^ + 2ixix 



\x)i\x + y + X2) dxdy 

= \Bal3ixi,X2))-i^2 ■ 



(22) 



A. Mathematics 



We will now proceed in calculating the entanglement 
of the first class of states, \Bai3{xi,X2))i2- We have al- 
ready shown that the entanglement should be equal when 
calculated using either partial density matrix (from the 
fact that a Schmidt basis exists) but we will show this 
explicitly by showing that the eigenvalue equations they 
give can both be transformed into the same eigenvalue 
equation. The two eigenvalue equations are 



exp 



1 



-' + ^) (-'+-") 



+2ix,ix - x') + ^) <j>^\x')dx' =AfVr^(^) (24) 



exp 



2(a2/?2 + 1) 



y^ + y'"^ 

0" 



yy 



c^f\y')dy' ^\fUf\x). (25) 



We can now recast the eigenfunctions of equation (p4|): 



Yi (^') ~^ Vi {x')exp{~2ixix') 



(26) 



thereby absorbing the complex exponential into the 
eigenf unctions. Next we perform the following change 
of variables in equation ( |24|) 



\/2x/3cr 



V2x'Pc 



(27) 



and in equation (25) 



X -^ x^j2{a'^fP + l)/3cr 



x' -^ V2(a2/?2 + \)x'iiG. 
(28) 



These changes of variables are the same for both dashed 
and undashed variables so this change can be absorbed 
into the eigenfunctions. All these changes leave the set 
of eigenvalues unchanged and give the same eigenvalue 
equation: 



exp 



(l + 2a^0^) (x? + a;'2) + 2xx'\ 0(a;') dx' 



K{x,x') 



Mix) (29) 



where K{x,x') is the Kernel of our integral eigenvalue 
equation. Its trace is 



K{x, x) dx 



2a(5 



(30) 



The set of eigenfunctions are of course independent of any 
normalisation constant but we will find that this constant 
is our primary check for the convergence of the numerical 
solution that follows in that the eigenvalues should sum 
to this constant as in equation (0). 

More importantly we notice that the eigenvalues are 
independent of the scale length, a^ and the value of xi 
or X2 ■ It is only dependent on the product of widths of 
the original Gaussian distributions with respect to a. 

The two particle cat state has partial density operator 
(aside from normalisation) 



1 1 

j=0 fc=o 

\x)^ {x'\ dxdx' , 



r 1 1 

ryy^ u ^*g-(:r-(-i)^d)'-(^'-(-i)'''i)'+2<i=5,fc 

^ A — n ;, — n 

(31) 



which is not easily transformed into a simpler form. Its 
trace is 



Tr{pi) = ^l-\e 



2d-' 



+ {AoAl + A*Ai)e 



-2d' 



(32) 



This, again, will be our primary check for convergence of 
the numerical model that follows. 



B. Numerical procedure for partially correlated 

states 



We cannot solve many of the integral eigenvalue equa- 
tions we encounter so we must use some numerical ap- 
proximation H. The most direct approach is to solve a 
discrete eigenvalue equation by approximating the inte- 
gral by the rectangle rule. Our eigenvalue equation has 
infinite limits so there must also be a cut-off point in the 
limits beyond which we do not approximate the integral. 

First, therefore, we discretize the eigenvalue equation 
(E9) into 2n + 1 parts (i = —n . . . , 0, . . . n) each of width 
5 covering the range —w < x < w where w — nS. Our 
eigenvalue equation then becomes 



^ X! Pp'i^P " ^'^9' 



(33) 



p— — 71 



where the indices now denote particular matrix and vec- 
tor elements (rather than particular eigenvectors) and for 
the Bell states 

Ppg = exp [(-(1 + 2«2^2)(p2 ^ ^2) ^ 2to) S^] . (34) 
Our entanglement measure is then 



^<-«-E(2^'-(^ 



(35) 



where the outer sum is over the set of eigenvalues and the 
sums over s and t are to normalise the set of eigenvalues. 

As mentioned above the primary check for the conver- 

tce of this solution is to compare ^ A with equation 
). We have two independent parameters that we may 
vary (at each value of a and /3) out of the three related 
parameters n, 6 and w. In practice for most values of a 
and f3, 2n + l — 201 and w around 10 standard deviations 
from the mean were sufficient for this sum to be equal to 
the trace accurate to 5 significant figures, this accuracy 
becoming greater with increased n and decreased d. 

We can see that what we are effectively doing in this 
procedure is sampling the spectrum of eigenvalues of the 
density operator over discrete ranges of width S by mod- 
elling the continuous system as a discrete 2n -1- 1 level 
system and taking the limit n — > oo. With the above 
convergence of the model we can be confident that both 
an adequate range of the integral has been sampled and 
that it has been sampled to an adequate precision. 

We used numerical procedures for eigenvalue problems 
from the NAG library to solve this problem for varying 
values of a and /3. The results are shown in figure [y. 
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FIG. 1. Entanglement of a partially correlated state in 
terms of a and /3, the widths of the Gaussians from which 
they are formed. 

As we expect the entanglement is increased when the 
parameters a and /3 are reduced, becoming infinite as 
these parameters approach zero, as will be shown in the 
next section. At this point it becomes harder to see con- 
vergence in the numerical procedure. 



C. Analytical results for partially correlated states 

We can now be even more general than this: we see 
that the form of the eigenvalue equation ( p9| ) is 

exp [-(1 + P){x^ + a;'2) + 2xx'] (j){x') dx' 

= \(i>{x) (36) 

with the entanglement being decreasing with increase in 
the parameter P. In fact, a state of the form 

exp [—ax — by + 2cxy + dx + ey) \x) \y) dxdy (37) 

has a parameter 

(38) 



, ab 
^ = 2(^-1 



independent of d and e, which can be shown by making 
appropriate linear changes of variables x — + Ax + B and 
x' — > Ax' + B in the eigenvalue equation formed by the 
density operators of the state ( |37| ) . 

For a two mode squeezed state with squeezing param- 
eter r 

^{x, y) = exp (^-i {e'^ix + yf + e-^^{x - yf)^ (39) 



P = 2cosech'^{2r). 



(40) 



Such a state can, of course, be written analytically in the 
number basis [Bj 



_. oo 

\^)i2 = —r7-:Y.i'anh{r)r\n), 



cosh(r) ^-^ 

^ ' n=0 



(41) 



which is the Schmidt basis for this state and from this 
the entanglement can be calculated as 

E = cosh {r)log2{cosh (r)) — sink (r) log2{sinh (r)). 

(42) 

It is reassuring to note that the above simulation can 
reproduce this analytical result accurate to 6 signifi- 
cant figures. More importantly we now have an ana- 
lytical result for the entanglement of the partially corre- 
lated states which is equation (B3) with the substitution 
2r — arcsinh[l/a(3), which follows directly from equa- 
tion M). 



D. Numerical procedure for two-mode cat state 

Now we move on to the entanglement of the cat states. 
The eigenvalue equation for these cannot be rewritten 
in terms of the parameter, P, so we proceed directly by 
discretizing the density operator of equation ( |3l] ) and 
making it the kernel of equation ( p3[ ) . Results for the en- 
tanglement of the cat states are shown in figure g against 
the parameters d and |v4op. 




FIG. 2. Entanglement of the cat states against (half) the 
distance between Gaussians, d, and the coefficient Aq. The 
entanglement is greatest for high values of d where the Gaus- 
sians become orthogonal, and for \Ao\'^ = 0.5 as with discrete 
entanglement. 



Note that the entanglement for any given value of d 
is maximum when |AoP = 1/2 and that the entangle- 
ment increases with d for given values of Aq, approach- 
ing the limit E = -\AQ\'^log2\Ao\'^ - \Ai\'^log2\Ai\^ as 
d —f oo, where the separated Gaussians become orthog- 
onal. Note also that only two eigenvalues dominated the 
contribution to the entropy, although we do not expect 
these states to be written in some basis as a two level 
system. These states, therefore, behave very much like 
discrete 2-level entangled systems. 



IV. PURIFICATION 



Alice 



Bell state projection' 



Bob 



initial 



final 



FIG. 3. The entanglement swapping procedure. Bob, hold- 
ing a copy of the entangled state shared by himself and Alice, 
performs a Bell state measurement on a particle from each 
pair and, for certain measurement outcomes, the entangle- 
ment of the final shared pair is higher than that of the initial 
shared pair. 



A. Imperfectly correlated states 



Now that we have a reliable method of calculating the 
entanglement we will move on to attempt entanglement 
purification in the two classes of states. The first point 
to note is that purification is always possible for such 
states as we can just make a measurement, one of the 
results of which is a projection onto the two levels of 
the Schmidt basis with the largest Schmidt coefficients. 
After this we can then perform discrete purification to 
obtain highly entangled two level states which may have 
higher entanglement that our original continuous state. 
However we wish to work more in the spirit of continu- 
ous systems using continuous operations and producing 
continuous entangled states. 

Again we generalize a purification procedure from dis- 
crete systems. The one we have chosen is purification by 
entanglement swapping [|0|. In the discrete case Alice 
and Bob share an entangled state and Bob holds another 
copy of the same entangled state as in the upper part of 
figure H. Bob performs entanglement swapping by mak- 
ing a Bell state projection on one particle from each pair 
of entangled particles. This leaves the remaining two 
particles (one on Bob's side and one on Alice's) in an en- 
tangled state. If the original entangled state were maxi- 
mally entangled then the final pair will also be maximally 
entangled. This is the standard entanglement swapping 
procedure |1^ . If the original state was less than maxi- 
mally entangled then the entanglement of the final pair 
will, for certain measurement outcomes, be less entangled 
than the original states, but for the remaining outcomes, 
will be maximally entangled. This procedure is, there- 
fore, probabilistic as there is dependence on the outcome 
of these measurements. Entanglement swapping has also 
been recently achieved independently in continuous vari- 
able systems [pSi. 



We will attempt the analogous procedure here except 
that the Bell state measurement will be replaced by a 
reverse entangling operation £^ and projective infinite 
resolution measurements on the two particles. 

The procedure for the partially correlated Bell state is 

2 {a\,{b\ (4 (|i?„„(0,c)),2 |i?^^(0,c))3j) 

ea^P _9^_9 , ^9^ ( -^^9{(^, P) - y^giP^ ol) + 2xy 

^)i\y)^ (43) 



_cr2(a2+/52) 

+2b{y -x)- 2ia(ya^ + x0^)\ 



where 



g{a,l3) 



a^l3' 



1. 



(44) 



Has the entanglement increased? First notice that from 
equation ( pq) the entanglement does not depend on a 
or b and so is not probabilistic. This indicates that the 
entanglement cannot have increased otherwise we would 
have deterministic purification. Indeed this is true as the 
parameter for this state is Pswap = 2[[a'^ + 0^ fi"^ + \){I3'^ + 
0^ a^ -\- \) — 1] whereas before the swapping process the 
parameter was Pq = 20* or 2/3^. But Pswap > -Fb in ei- 
ther case and the entanglement is strictly decreasing with 
increase in P so the swapped pair has less entanglement. 
Of course, our final projections in this method were 
onto the unphysical states \a) and \b) but further cal- 
culations indicate that with finite width projections the 
parameter P still increases. Such calculations involve 
Qth (jggj-ge polynomials in the width parameters {a and 
/3 etc.) so proving that P increases for all values of these 
parameters is difficult and we have not been able to do 
so analytically. However, graphical results indicate that 
this is true. 



B. Two- mode cat states 



We now attempt a similar method with the two-mode 
cat states, but setting the scale length a = I and making 
finite resolution measurements: 



/14 



2 (G^(a)l 3 {G^rn (4 (|C(rf))i2 |C(d))34) 



-(x-(-iydf-(y+(-i)''df 
^ g(dh/i(M)((-l)^+(l+M')(-l)'')+2d^5jfc) 

^ g(»ad/l(M)((l+M')(-lV-(-l)'=)) la;) \y 



dxdy (45) 




a/it 



FIG. 4. Entanglement of swapped cat states with ^ = 0. 
Above the level of the plane purification has been achieved. 



where 



hi^,) = 



2 + 2/^2 + ^4 ■ 



(46) 



Writing this state out in full in the high precision mea- 
surement limit, /i = 



w 



14 



AoAn e-2d'+2'^^ ^-{x-df-iy+df 



g2md ^~{x-d) -{y- 



-AqAi 

-AiAo e-2»ad ^-{x+df-iy+df 

-AiAi e-^d'-adh ^^{x+dy'-iy-d)^ 



\x), |j/)4 dxdy (47) 



and looking at the particular case where the probabilistic 
measured values are a = 6 = we can see purification 
for high values of d as the middle two terms now domi- 
nate and have equal coefficients. As d — > oo they become 
maximally entangled. This is again very much like the 
action of discrete entanglement under purification proce- 
dures: the coefficients of the states have changed not the 
states themselves. 

For the results of figure and we have chosen the 
values ^0 = VO^ and d — 1.0. They show the entangle- 
ment of the resulting state (Esh for a range of values of a 
and b with /i = and 0.5 respectively. 




aAr 



FIG. 5. Entanglement of swapped cat states with n — 0.5. 
Again purification has been achieved above the level of the 
plane, but with the inaccuracy in the measurement part of 
the entanglement swapping process the amount of purifica- 
tion is reduced. 

The horizontal planes are at the level of entanglement 
of either cat state before the purification procedure is 
performed, that is, above this plane purification has been 
achieved. Notice that purification is achieved above the 
level of entanglement of the initial cat state of equation 
( p3| ) with parameters Aq ~ VO.S and d = 1.0 for which 
^ = 0.881. 



V. CONCLUSIONS 

We have shown that the usual measure of entanglement 
for pure bipartite discrete systems is easily generalized to 
continuous systems with knowledge of the mathematical 
area of integral eigenvalue equations, in particular that a 
Schmidt basis exists provided the state satisfies certain 
reasonable conditions of integrability. The calculation of 
the entanglement, however, is often difficult so a numeri- 
cal procedure was necessary. The results of the numerical 
procedure corresponded very well with those of analyti- 
cal results where these were available and enabled us to 
test a purification procedure, the entanglement swapping 
procedure, for the two classes of states we presented. 

Curiously this procedure only succeeded for one of the 
classes of states, the two-mode cat states whose charac- 
teristics are very much like those of discrete systems, al- 
though there does not appear to be a discrete and finite 
basis in which the states can be written. The purification 
procedure also acted in a similar manner to the analogous 
discrete procedure. 

For the other class of states, the partially correlated 
states, no continuous procedure could be found that in- 
creased the entanglement in the state, indeed, no proce- 
dure was found where the entanglement of the final state 
was in any way dependent on the measurement results 
during the procedure, reassuring us that no purification 
would be possible. However, we plan to study this prob- 
lem more carefully. 

What exactly the key difference is between these two 
types of states which allows purification in one class but 
not the other is still unclear. There are obvious cor- 
respondences between the form of entanglement in the 
two mode cat states and discrete systems and it would 
be interesting to find a condition for continuous variable 
purification, as it has been attempted here, which a state 
undergoing purification must obey. The fact that purifi- 
cation has been demonstrated, however, in continuous 
systems is an interesting result. 



ACKNOWLEDGEMENTS 

This work is supported by the United Kingdom En- 
gineering and Physical Sciences Research Council (EP- 
SRC), the Inlaks Foundation, The Leverhulme Trust, 



the EU TMR-networks ERB 4061PL95-1412 and ERB 
FMRXCT96-0066 and the European Science Foundation. 



[1] S. L. Braunstein, H. J. Kimble, Phys. Rev. Lett. 80, 869 

(1998). 
[2] G. J. Milburn, S. L. Braunstein, Quantum tel eporta 



tion with sq ueezed vacuum states, LANL eprint: |quant- 
ph/98120l4 



[3] S.J. van Enk, A discrete formulatio n of teleportation of 

continuous variables, LANL eprint: quant-ph/9905081. 
[4] T. C. Ralph, P. K. Lam, Phys. Rev. Lett. 81, 5668 (1998). 
[5] A. Furusawa, J. L. S0rensen, S. L. Braunstein, C. A. 
Fuchs, H. J. Kimble, E. S. Polzik, Science, 282 23 Octo- 
ber, 706 (1998). 
[6] U. Leonhardt, Measuring the Quantum State of Light, 

Cambridge University Press 
[7] C. H. Bennett, H. J. Herbert, S. Popescu, B. Schumacher, 
Phys. Rev. A 53, 2046 (1996); S. Popescu, D. Rohrlich, 
Phys. Rev. A 56, R3319 (1997). 
[8] D.A. Kofke, A. J. Post, J. Chem. Phys. 98, 4853 (1993). 
[9] D. Deutsch, A. Ekert, R. Jozsa, C. Macchiavello, S. 
Popescu, A. Sanpera, Phys. Rev. Lett. 77, 2818 (1996); 
C. H. Bennett, D. P. DiVincenzo, J. A. Smolin, W. K. 
Wootters, Phys. Rev. A 54, 3824 (1996). 
[10] S. Bose, V. Vedral, P. L. Knight, Phys. Rev. A 60, July 

(1999). 
[11] S. Braunstein, Phys. Rev. Lett. 80, 4084 (1998). 
[12] F. G. Tricomi, Integral Equations, Wiley, New York 

(1967). 
[13] A. Wehrl, Rev. Mod. Phys. 50, 221 (1978). 
[14] L. Vaidman, Phys, Rev. A 49, 1473 (1994). 
[15] B. C. Sanders, Phys. Rev. A 45, 6811 (1992). 
[16] C. C. Gerry, Phys. Rev. A 55, 2478 (1997). 
[17] M. Zukowski, A. Zeilinger, M. A. Home, A. K. Ek- 
ert, Phys. Rev. Lett. 71, 4287 (1993); J -W. Pan, D. 
Bouwmeester, H. Weinfurter, A. Zeilinger, Phys. Rev. 
Lett. 80, 3891 (1998); S. Bose, V. Vedral, P. L. Knight, 
Phys. Rev. A 57, 822 (1998). 
[18] R. E. S. Polkinghorne, T. C. Ralph, Entanglemen t swap 



ping using continuous variables, LANL eprint: quant ^ 
ph 9906066]; P. van Loock, S. Braunstein, Unconditional 



entangl ement swapping fo r continuous variables, LANL 
eprint: quant-ph 9906075. 



